Electric Control of Spin Currents and Spin- Wave Logic 
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Spin waves in insulating magnets are ideal carriers for spin currents with low energy dissipation. 
An electric field can modify the dispersion of spin waves, by directly affecting, via spin-orbit coupling, 
the electrons that mediate the interaction between magnetic ions. Our microscopic calculations 
based on the super-exchange model indicate that this effect of the electric field is sufficiently large 
to be used to effectively control spin currents. We apply these findings to the design of a spin- 
wave interferometric device, which acts as a logic inverter and can be used as a building block for 
room-temperature, low-dissipation logic circuits. 



One of the major challenges of contemporary electron- 
ics is to reduce dissipation as the size of devices shrinks to 
the nanometric scale. In this context, spin- wave spintron- 
ics, so called magnonics, with insulating magnets offers 
interesting possibilities [l|, |2[ . While in metals and semi- 
conductors the spin current is carried by mobile conduc- 
tion electrons/holes, which inevitably dissipate energy as 
they move, in a magnetic insulator, such as Y^Fe^0i2 
(YIG), the spin current is carried by a collective motion 
of magnetic moments - a spin wave - with no charge dis- 
placed. The spin current propagating in these insulating 
material is thus totally free of energy dissipation from 
Joule heating, and almost free of dissipation from other 
sources (e.g. electron-magnon scattering) : the coherence 
length can be as large as several centimeters [Ij. For 
these reasons, magnetic insulators have attracted con- 
siderable attention in recent theoretical 0^ and experi- 
mental 0, 0] work. For example, Kajiwara et al. , have 
demonstrated injection and extraction of spin waves into 
and out of a YIG wave guide [§]. Kostylev et al. @, 
have designed an ingenious scheme of spin- wave logic, 
based on the interference between spin waves traveling 
along different arms of a Mach-Zehnder interferometer 
(a schematic illustration of a Mach-Zehnder spin wave 
interferometer is shown in Fig. [T]). 

A crucial element of magnonics [l| is the phase shifter 
- a device that changes the phase of propagating spin 
waves. Several mechanisms have been proposed in the 
past to implement controlled phase shifts on spin waves. 
The simplest and most direct, is the application of a mag- 
netic field, which shifts the dispersion jlO|. thus changing 
the wave vector at constant frequency [9|. More sophis- 
ticated mechanisms exploited the Berry phase accumu- 
lated by spin waves that propagate on a non-collinear 
magnetic texture 0, El- In a parallel development, Cao 
et al. 



Figure f. A Mach-Zehnder spin- wave interferometer in the 
presence of radial E field. A weak magnetic field is applied 
perpendicular to the ring plane, tilting the equilibrium mag- 
netization away from the ring but still in the tangential plane 
to the ring. Oq denotes the orientation of the equilibrium 
magnetization. 



current) by means of an electric field. We will show that 
the electric field-induced AC phase has important im- 
plications for spin wave interferometry, since the effect 
is much larger than that initially predicted in Ref. |12 . 
Our analysis starts at the microscopic level, with a very 
simple super-exchange model (TgI for the magnetic in- 
teraction between two neighboring magnetic ions (e.g. 
Fe'^"'" in YIG) in a magnetic insulator. The model is de- 
picted in Fig. [21 There are no itinerant electrons in an 
insulator, but the virtual hopping of electrons between 
the d-orbitals on the magnetic ions (Fe'^"'') and the p- 
orbitals on the the ligand (0^~) is sufficient to establish 
an antifcrro-magnctic interaction of the Heisenberg type. 



Hh = JSi ■ S2 



(1) 



12| studied the effect on spin waves of an electric 
field-induced Aharonov-Casher (AC) phase More 
recently, the influence of electric fields on spin waves 
has been studied both theoretically [l^ and experimen- 
tally [l^ and a strong shift of spin-wave dispersion in- 
duced by an electric field has been reported [1^ . 

In this Letter we directly tackle the problem of con- 
trolling the phase of a spin wave (and hence the spin 



which is responsible for the occurrence of magnetic order 
in the material. Notice that in this model the physical d 
and p orbitals are replaced by doubly degenerate orbitals, 
which are eigenstates of the z-component of the electron 
spin. Spin-orbit coupling effects are completely neglected 
up to this point. 

We now augment the usual super-exchange model by 
the inclusion of a spin-orbit (SO) interaction of the form 



A 

Hso = — r(P X eE) • cr , 

n 



(2) 



where A is a characteristic length scale that controls the 
strength of the SO interaction, e is the elementary charge 
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Figure 2. Super-exchange modeh two half-filled magnetic ions 
connected by an oxygen ligand. 



and cr is the Pauli matrix. For electrons in vacuum A is 
the Compton wavelength Ac = , but we will see below 
that, in any realistic model of magnetic insulators, the 
value of A is orders of magnitudes larger (of the order 
of A). Although there are no itinerant electrons carry- 
ing a finite average momentum p, we will show below 
that the effect of the SO interaction on the phases of the 
virtual hopping translates, at the macroscopic level, into 
the appearance of a Dzyaloshinskii-Moriya (DM) inter- 
action [17[ between the magnetic moments of the ions: 



Hdm = D • (Si X S2) 



where the vector D is given by 



D 



-J— — E X ei2 ; 



(3) 



(4) 



with a being the distance between the magnetic ions, 612 
the unit vector in the direction connecting the ions, and 
Esn = > (in vacuum A Ac and Eso ~ 



2mA2 



m is the bare electron mass). As a result of this interac- 
tion, the spin waves - obtained by solving the appropriate 
ferromagnetic Hamiltonian for an infinite chain of identi- 
cal blocks of magnetic ions (see discussion below) - have 



their wave vector shifted by q 



D Sc 

J a 



ei2, where So is 



the direction of the equilibrium magnetization. This, in 
principle, gives us a way to control the phase of the spin 
wave by an electric field. In practice, the feasibility of the 
proposal depends critically on the strength of the spin- 
orbit coupling A^. If we used the value of A in vacuum 
(as was done in Ref. [12), the effect would be extremely 
weak, and probably unobservable with realistic electric 
fields. 

The reason why A turns out to be much larger than 
Ac is that the physical d-orbitals in the magnetic insu- 
lator (as opposed to the model orbitals we have dealt 
with so far) have strong intrinsic spin-orbit coupling L • S 
built in. For example, in the model adopted by Katsura 
et al. 0, the doubly degenerate orbitals of the super- 
exchange model are actually spin-orbit-entangled states 
of the form \a) = {\dxy t) + \dyz i) + A^zx i))/V3 and 
\b) = {\dxy i) - \dyz ■\)+i\dzx t))/A/3- In order to deter- 
mine the value of A (or, better, Eso) in our model, we 
observe that the DM interaction Eqs. p|4p . and indeed 
the spin-orbit interaction ^ itself, can be viewed as the 
interaction of the electric field with an effective electric 



dipole, i.e. Hom = E • P, where the electric dipole is 
eo 

P = -./i^ei2 X (Si X S2) . (5) 

The effective electric dipole arises from the hybridiza- 
tion of orbitals centered at different atoms. For example, 
in the model of Ref. Q the exact single-particle eigen- 
states arc combinations of d orbitals on the magnetic 
ions and p orbitals on the oxygen due to the hopping of 
the electron: these states carry an electric dipole moment 
~ ~^M^^i2 ^ (^1 ^ ^2) , where t is the hopping coeffi- 
cient and / = if Zo'^^zjf (^ + ^)-^aB with as being 
the Bohr radius and Zq {Zm) being the atomic number 
of O (M). Comparing this to our Eq. ([5]) we arrive at an 
unambiguous identification of Eso (and hence A) within 
our model: 



9ta 

Eso = ^. 



(6) 



Taking YIG as an example 3 3, with t = 0.8 eV, and 
/ = 0.61a, we get Eso = 3.0 eV and A = 1.13 A. This 
is indeed several orders of magnitude larger than Ac and 
opens the way to practical schemes of electric control of 
the phase of spin waves. 

In the remaining part of this paper we supply more the- 
oretical detail on the calculations supporting the above 
analysis, then work out the dispersion of spin waves in 
the presence of the electric field, and apply the results of 
these calculations to the design of a spin- wave interferom- 
eter (see FigH]), to be used as building block of spin- wave 
logic circuits. 

Microscopic analysis - The insulator we are interested 
in is YIG, whose magnetic order mainly arises from the 
super-exchange interaction between Fe^~^ in octahedral 
(a) sites and tetrahedral (d) sites. The fact that the 
numbers of (a) and (d) sites per unit cell are different 
makes YIG a ferrimagnet. However, the long wave-length 
spin waves, whose energy is less than about 40 K, can 
be understood with an effective ferromagnetic exchange 
coupling between "block spins" Si, one per unit cell [l9|. 
The question is how an electric field affects this block 
ferromagnet. We start from the super-exchange model 
shown in Fig. [5J which can be described by the following 
Hamiltonian 

H super = Ho + Ht + Hu , 



Ht = -t^{c{^CQa + 4<tC0(t + h.c.) , 



(7) 



where c (c^) is the creation (annihilation) operator of 
ligand electrons, which can hop forth and back only be- 
tween oxygen ligand and the metal ions, ei and ep are 
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the orbital energies of a metal ion and the oxygen lig- 
and, respectively. The large repulsion energy [/ (~ 8 
eV) between two electrons on the same metal ion allows 
for a maximum occupancy of two electrons per ion (the 
repulsion between the electrons in the oxygen ligand is 
negligible in comparison). 

The fact that t (~ 0.8 eV) is much smaller that U al- 
lows us to use perturbation theory. Keeping up to the 
fourth order of t yields the effective interaction between 
the spins on the magnetic ions: 




where V = ei — eo + U corresponds to the energy dif- 
ference between the p and \P)j orbitals in the paper by 
Katsura et al. . Setting J = + ^ ^ and drop- 
ping the constant term, we obtain the Heisenberg inter- 
action Eq.(IT]). A positive J implies that the interaction 
between neighboring magnetic ions is antiferromagnetic. 
However, this antiferromagnetic interaction gives rise to a 
ferromagnetic interaction between "block spins" in YIG, 
due to the unequal magnitudes of the anti-parallel mag- 
netic moments in each block. 

Let us now include the spin-orbit interaction Hso from 
Eq. It is easy to see that the inclusion of this inter- 
action is equivalent to the inclusion of a spin-dependent 

,2 

vector potential A = ■^^^"-^ ^ ''^hicli in turn modifies 
the hopping term Ht by a spin-dependent phase factor, 
that is 

Ht = ^(cLco.e-*"'^ + cLco.e*"" + h.c.) , (9) 

fT 

where a = -^f^, provided that the external electric field, 
the motion of the electron, and the electron spin (ct) are 
perpendicular to each other. Notice that the phase a is 
proportional to the distance between neighboring sites 
and independent of the direction of the localized mo- 
ments: one can therefore switch to the "block spins" 
description by simply reinterpreting a as the distance 
between neighboring blocks. The resulting spin Hamil- 
tonian takes the form 

H=-.r J2 S^S] + ^{e'^"^^S+S- +e-'^"^^S-S+) 
^ -J' {(S. • S,) +sin2a,,(S, x S^-).} , (10) 

<ij> 

where atj = 2a(i — j), —J' is the effective exchange 
coupling for the spin blocks and z is in the direction 
perpendicular to E and e^j. In addition to the nor- 
mal Heisenberg term we now also have a DM term [l7j . 
whose strength is linear in E. An electric-field induced 
anisotropy is also present, but is an effect of order i?^ 
and has therefore been neglected for weak electric field. 

In spite of the presence of the noncollinear DM term, 
the ferromagnetic configuration is still the ground state 



of (|T0)) . To show this, we make Si = Sq + 5Si, where 
dSi is small deviation perpendicular to Sq- Then the 
variation of the DM term up to the second order of SSi 
is 

<i,j> 

where D^^- is defined by Eq.Q with ei2 replaced by e^. 
Since = —eji, we see that D^j — 0, which means 
SHdm = up to the first order of 6Si. Hence, the ground 
state is still ferromagnetic. However, the DM term will 
definitely modify the spin-wave frequency, which involves 
a correction to the ground state energy at the second 
order in (5Sj. Further, we can clearly see that it is only 
the component of the D parallel to ^S^ x SSj (i.e. to the 
direction of the equilibrium magnetization) that plays a 
role in the modification. 

Spin wave dispersion - We now proceed to solve the 
dispersion of the spin waves in the presence of the DM 
interaction derived above. We consider the ring geometry 
illustrated in Fig[T] the electric field perpendicular to the 
ring produces a DM vector D directed along the z-axis. 
This will affect the dispersion of spin waves if and only if 
the equilibrium magnetization has a non-vanishing com- 
ponent along the z axis. In a flat ring, such as the one 
shown in FigdJ the shape anisotropy [2lJ —K{S ■ e)^/S'^ 
where e is the unit vector along the ring - outweighs other 
forms of anisotropy, causing the equilibrium magnetiza- 
tion to lie along the ring, in which case the electric field 
has no influence. As a result, a magnetic field along the 
z axis (Zeeman coupling g^BBSz) is necessary for us to 
observe the impact of the DM term on spin waves prop- 
agating in the ring. Now, however, the orientation of the 
equilibrium magnetization is no longer constant in abso- 
lute space (even though it is constant relative to the ring). 
This causes an additional geometric phase (ag = ■^) to 
appear, as shown in Ref. 0, where R is the radius of the 
ring. Putting everything together, i.e., DM interaction, 
geometric phase, Zeeman coupling and shape anisotropy, 
we arrive at the following equation of motion: 

= S, X [J'(S,+i + S,_i) + 2^5fe, - ffMB-B.e,] 

-i?,5f(S,+i -S,_i). (12) 

The large magnitude of the "block spin" of YIG (S=14.3) 
allows us to use the semiclassical spin-wave approach to 
get the dispersion relation: 

J' s I '. 

cj = -j^{a^\J (fc2 + K2)(fc2 ^ K^sin^ 6*0) + 2afcacos6'o} , 

(13) 

where a = Aa — ag, k = \J ji^i^'i 1 s-nd cos Oq = ^^^"^ , 
which is determined by minimizing the total Hamiltonian 
in the limit of 3> 

As shown in Fig. [31 one can tune the dispersion by 
adjusting the electric and the magnetic fields. Just as a 
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Figure 3. (a) Dispersion of spin waves in the ferromag- 
netic ring in Fig[T] taking into account the geometric phase 
and the phase induced by the electric field. Parameters we 
use [H: J' = 1.18 x 10"" eV, K = 1.53 x 10"* eV, S= 14.3, 
a = 12.4 A, ro = 50 nm, R = 100 nm, B = 0.05 T. (b) Trans- 
mission probabihty of a spin wave in the ring interferometer as 
a function of input vohage and magnetic field at = 43 GHz. 



magnetic field shifts the spin wave dispersion vertically 
by increasing or decreasing the frequency at fixed fc, the 
electric field shifts the dispersion horizontally by increas- 
ing or decreasing the wave vector at fixed frequency. 

Spin wave interferometer - Now we are ready to de- 
sign our spin-wave interfcrometric device. An insulating 
ring encircles a metal electrode to which a voltage Vin 
can be applied. The radial electric field acting upon the 
electrons in the ring is pi"X'"/ps . 

In Fig. [3] (b) we plot the transmission of a spin wave 
sent through this Mach-Zender interferometer, as a func- 
tion of Vin and B. The effect of B is to change the equilib- 
rium orientation of the magnetization. The white regions 
in the figure are regions of constructive interference, sep- 
arated by regions of destructive interference. We see that 
very modest changes of potentials and magnetic fields, of 
the order of 1 V and 0.01 T respectively, switch the re- 
sponse of the interferometer from high to low. It is then 
clear how the device can be used as a logic inverter: the 
logic input being the voltage on the central electrode, 
and the logic output the intensity of the spin wave, as 
measured by an inductive coupler. Advantages of this 
design are that it would operate at room temperature 
and GHz frequencies, with very little dissipation, and 
can be made small by using exchange spin waves - the 
only type we are really considering here, since magneto- 
static spin waves have much longer wavelengths and are 
hardly affected by the AC phase. Once a logic inverter is 
available, we can follow Kostylev et al. in construct- 
ing more complicated architectures, which implement the 
NAND, the NOR, and all of classical logic. 

It is worth noting that a traveling spin wave is itself 
a source of electric field: E ~ /iov x (5M where v is the 
velocity of the spin wave, (5M is the amplitude of the 
magnetization oscillation and /xq = 47r x 10~^N/A^ is 
the vacuum permeability. In our device, the resulting 
electric field is of the order of 1 V/m, which is negligible 
in comparison to the control field E ~ 10^ V/m. Even 
smaller is the electric field associated with mesoscopic 



equilibrium spin currents in the ring (E ~ 10"^ V/m). Q 
In conclusion, we have proposed an energy-efficient 
way to control the spin current propagating in an insulat- 
ing magnet by means of an electric field. This possibility 
arises from the strong coupling that exists between the 
electric field and the spins of the electrons that medi- 
ate the interaction between magnetic ions. The strength 
of this coupling has been theoretically estimated from 
microscopic parameters, such as electron hopping coeffi- 
cient, distance between neighboring magnetic sites, etc. 
as shown in Eq.®. Or, it could be indirectly determined 
from measurements of physical effects that are sensitive 
to it, e.g. the spin wave spin Hall effect proposed by 
Meier and Loss [5[. Finally, we have applied our the- 
ory to an insulating magnetic ring inteferometer, which 
can be used to implement a voltage-controlled spin- wave- 
based NOT gate. 
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